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Abstract 

We use Ikehata's enclosure method to reconstruct penetrable unknown in- 
clusions in a plane elastic body in time-harmonic waves. Complex geometrical 
optics solutions with complex polynomial phases are adopted as the probing 
utility. In a situation similar to ours, due to the presence of a zeroth order 
term in the equation, some technical assumptions need to be assumed in early 
researches. In a recent work of Sini and Yoshida, they succeeded in abandon- 
ing these assumptions by using a different idea to obtain a crucial estimate. 
In particular the boundaries of the inclusions need only to be Lipschitz. In 
this work we apply the same idea to our model. It's interesting that, with 
more careful treatment, we find the boundaries of the inclusions can in fact 
be assumed to be only continuous. 

Keywords: enclosure method, reconstruction, complex geometrical optics 
solutions, time-harmonic elastic waves. 



1 Introduction 

In this paper we consider the inverse problem of reconstructing penetrable unknown 
inclusions in a plane elastic body by boundary measurements. In [T7] and [19], the 
same problem is considered in the context of elastostatics. In the present work we 
shall consider the situation when time-harmonic waves are applied. The mathemat- 
ical model is described in the following. 



1.1 Mathematical model 

Let C be a bounded domain (open connected set) occupied by our object, which 
consists of an elastic body as background and some unknown inclusions therein. For 
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simplicity we assume Q has C°° boundary. The background elastic body will be 
assumed to be homogeneous and isotropic with Lame constants denoted by Aq and 
/iq. Denote the region of unknown inclusions hy D. D is an open subset of with 
D (Z ^. The inclusions are also assumed to be isotropic but may be inhomogeneous. 
Denote the differences between the Lame coefficients of the inclusions and the back- 
ground by A/) and fi£,, which are assumed to be in L°°{Q), with \d = A^d = on 
fl\D. So the Lame coefficients A and /i of the whole object on fl are given by 

A = Ao + Ad and // = //o + A*d- 

For simplicity we also assume our object has unit density. Now, consider we send 
a time-harmonic elastic wave with time dependence e*'^* into Q. By singling out the 
space part wc have the displacement field u, which is a two-component vector- valued 
function, satisfying 

V • ((7(u)) + A;^u = mQ. (1.1) 

Here, for any displacement field v (which we will assumed to be a column vector), 
a(v) is the corresponding stress tensor, which is represented by a 2 x 2 matrix: 

a{v) = A(V • v)/2 + 2//e(v), 

where I2 is the 2x2 identity matrix and e(v) = |(Vv -|- (Vv)^) denotes the in- 
finitesimal strain tensor. Note that for v = (^1,^2)"^; Vv denotes the 2x2 matrix 
whose j'-th row is Vvj for j = 1,2. And for a 2 x 2 matrix function A, V ■ A denotes 
the column vector whose j-th component is the divergence of the j-th row of A for 
J = 1,2. 

For D — that is for the case with no inclusion, the corresponding displacement 
field will usually be denoted by Uq, which satisfies 

V ■ ((7o(uo)) + A;^uo = in Q, (1.2) 

where 

ao(v) = Ao(V ■ v)/2 + 2/ioe(v) 

for any displacement field v. Accordingly, we will use aoi'v) to denote cr(v) — cro(v), 
i.e. 

ctd(v) = Ac(V • v)/2 + 2/ioe(v). 
We assume Xq, Ho and A,/i satisfy the conditions 

Ao + 2/io > 0, /io > 0, and 

A + 2/X > 0, ^ > on O, ^ ' ^ 

which ensure respectively that —V • (Tq and —V • a are strongly elliptic operators. 

In particular the two operators both have at most countably many Dirichlet eigen- 
values. As a consequence, we can readily choose (and will choose) A; e M so that 
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is neither an eigenvalue of —V • ctq nor an eigenvalue of —V ■ cr. In this situation, 
the Dirichlet boundary value problems corresponding to fll.ip and (11.21) have unique 
solutions (see e.g. Ch.4 of [H]). Thus we can define the Dirichlet-to-Neumann maps 
Ad and A0, both from H^dnf to R-^Ony, by 

Az)f = a{u)iy\Qn and Ami = o-o(uo)i/|oo, (1.4) 

where u is the unit outer normal on dfl and u and uq solve respectively (II. ip 
and (11.21) with Dirichlet boundary data f . The goal is to determine the unknown 
inclusions from the knowledge of and A0. 

1.2 The method and improvement 

We will utilize the enclosure-type method to reconstruct the unknown inclusions. 
Such kind of methods are initiated by Ikehata and have been successfully applied 
to a various type of reconstruction problems, see for example [H El El El- In this 
method, complex geometrical optics (CGO) solutions usually play the important role 
of the probing utility. In his early works, Ikehata use the Calderon type harmonic 
function [T] e^ *^'^"'"*'^ ^ with u G S"""^ [6l[5]. It looks like one uses lines (planes) to 
enclose the obstacle (and hence the name). As a consequence a connected inclusion 
is required to be convex for a complete identification, and in general only its convex 
hull can be determined. One can refer to the survey paper [7j for detailed explanation 
and early development of this theory. In [15], [H] and [3], the authors utilize the 
complex spherical wave solutions and some concave parts of unknown inclusions can 
be determined. In Q^, the authors proposed a framework of constructing CGO 
solutions with general phases for some elliptic systems in two-dimension. In the 
same paper they then applied CGO solutions with complex polynomial phases to 
conductivity equations, and inclusions with more general shapes can be determined. 
This type of CGO solutions were later applied to other equations, for example [19] 
for static elastic systems and [131 fo'^ Helmholtz equations. In this work we will also 
apply CGO solutions with complex polynomial phases to our problem, of which the 
governing equations are the Helmholtz type elastic systems (11.11) and (II. 2p . 

A crucial point in our problem, as in [HI [TH [13], is the presence of the zeroth 
order term. Due to this, some technical assumptions are needed in early researches. 
In particular dD is assumed to be C^. However in the recent work [I6l of Sini and 
Yoshida, by using a different idea to obtain a crucial estimate, they succeeded in 
abandoning these technical assumptions, and in particular dD can be only Lipschitz. 
In this paper, we apply the same idea to our model. With more careful treatment, we 
find the boundaries of the inclusions can in fact be assumed to be only continuous. 
More detailed discussions are given in the remark after our main theorem. Theorem 

In the following we give a sketch of this paper as well as a rough idea of the 
whole process of the enclosure method. In section 2, we introduce a functional E 
on H2(^dQ)'^, which will be called the indicator functional in this paper. And then 
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we give an upper bound and a lower bound of E, which play central roles in the 
proof of the main theorem. In fact, we will construct a family f^/^ G H^ldQ)"^ as 
input data into E, and the limiting behavior of the output data, for various d, will 
indicate the location of dD. The construction of id,h is based on the construction 
of CGO solutions for (11.21) . which is given in section 3. By using the Helmholtz 
decomposition and the Vekua transform, this construction is much the same as in 
[13]. The main theorem concerning the limiting behavior of E on f^ /^, as well as 
discussions on the implication, the idea of proof and our improvement, are given in 
section 4. 



2 The indicator functional 

In this section we introduce the functional E on (dVl) defined by 



m 



[{Ad - A0)f] ■ ids, 



on 



where the Dirichlet-to-Neumann maps and A0 are defined in (11. 4p . E will be 

called the indicator functional (according to Ikehata's indicator function [5]), which 

plays a central role in the enclosure method. Intuitively, it measures, for a fixed 

Dirichlet boundary data, the difference between the tractions corresponding to the 

situations with and without D. 

Now let u and uq G H^{Qy satisfy (II. ip and (ll.2p respectively with the same 

boundary condition f G H^[dQ)'^, and let w = u — uq. The goal in this section 

is to prove Lemma (2.21 which gives a lower bound and an upper bound of E{i) in 

terms of uq and w. To this end, we first give two identities. Note that we use |y4| 

/ \ 1/2 

to denote ( ■ af^ j for a matrix A = (flij). 

Lemma 2.1. We have the following two identities: 



^|(AD + /iD)|V-uo|' + 2/i 



e(uo) - -(V ■ Uo)/2 



dx 



e(w) - -(V • w)/2 



dx 



(2.1) 



k"^ |w| dx; 



m 



J |(AD + /ii5) |V-u|V2/ifl 
+ ^|(Ao + /io)|V-w|' + 2/io 



e u 



(V ■ u)/2 



dx 



dx 



(2.2) 



k Iwl dx. 
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Lemma 2.2. Assume that the Lame coefficients Aq, /Xq one? A, jj, satisfy the strong 
convexity condition, that is 

Ao + /io, /^o > and X + fi, > 0, 
then we have the following upper bound and lower bound of E{i): 

Eii) < [ (Xd + f^D) \^ ■ uof dx 
Jd 

+ 2 / /ifl 



m > [ 

JL 



D 



dx + I k'^ \wf dx] 
n 



(Ad + /Ud)(Ao + /Uo) I |2 , 
|V ■ uqI dx 



D X + fl 



+ 2 



D 



e(uo) - -(V • Uo)/2 



dx — I k'^ Iwl^ dx. 



Proof. The upper bound of E{i) follows immediately from (12. ip (by omitting the 
second integral). On the other hand, from (12. 2p we have 



^(f)>y|^|(AD + M|V-u|2 + 2^ 



D 



+ J |(Ao + /io)|V- w|' + 2/io 

— / k'^\w\'^dx. 
Jn 



e(u)--(V-u)/2 



dx 



dx (2.3) 



And the lower bound follows from the following two identities, of which the verifi- 
cations are straightforward (by using w = u — uq). 



(i) 



(Xd + /Ud) |V ■ u|' + (Ao + /io) |V ■ w|' 

= fv^V-u-^^V-uoy + 

V V A + /i / 



(AD + /iD)(Ao + /io) 1^ |2 

^— V ■ uo 

A + 



2/x 



D 



e(u) - -(V . u)/2 



E 



+ 2/io 

2 



e(w) - -(V ■ w)/2 



e(uo) - -(V ■ uo)/2 



where 



(6.,):=e(u)--(V-u)J2 and m := e{uo) - -{V ■ Uo)l2. 



5 



□ 



For completeness we give the proof of Lemma 12.11 in the following. Before doing 
so, note that we have the following basic formulae: 



V ■ (ct(u)v) = (V ■ or(u)) ■ V + tr((T(u)Vv); 
tr(a(u)Vv) = tr((T(v)Vu). 



Here tr(-) is the trace of matrices. And 

tr(a(u)Vu) = (A + /i)|V ■ up + 2/i 



6(u)--(V-u)/2 



(2.4) 
(2.5) 



(2.6) 



These formulae are easy to check and we shall omit the proof. Also note that we 
have similar formulae with a replaced by do, ob-, etc. 
Now we give the proof of Lemma 12.11 

Proof of Lemma \2.1\ First note that /^^ A^jf • ids and /^^^ A0f • ids are real. In fact, 
by definition we have 



A/)f ■ ids = / ((t(u)i/) ■ uds 

an Jan 



cr(u)"^u) ■ vdx. 



an 



By divergence theorem and (12. 4p we then get 



ADf-fds= / (y ■ a{u)) ■ udx + / tr{a{u)Vu)dx 

Jn Jn 

which is real. Similarly J^^ A0f ■ ids is real. 

Since u and uq both equal f on dQ, similar to fl2.7p we have 

/ Aj;)i-ids= / —k^u-UQdx+ / tr{a{u)Vuo)dx; 
Jan Jn Jn 

/ Agf ■ f c?s = / —k'^Uo-udx+ / tr(cro(uo)Vu)(ia:. 
Jan Jn Jn 

Take complex conjugation of (12.81) and by (12.51) we get 

/ A£)i-ids= / —k^UQ-udx+ / tr(a(uQ)'Vu)dx. 
Jan Jn Jn 

Then subtract (12.91) from (12.101) we obtain 



(2.7) 



(2.8) 
(2.9) 

(2.10) 



E{i) = / tr{aDino)Vu)dx. 
'n 



(2.11) 
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On the other hand, 

/ /c^w • vfdx = I (A;^u — /c^uo) ■ wdx = — / V ■ (cr(u) — cro(uo)) ■ wdx. 
Jo, Jn Jn 

Note that w G Hq{Q)'^, thus integration by parts gives 

k"^ \wfdx= / tr [(a(u) -ao(uo))Vw](ia;. (2.12) 
Jn Jn 

Now, substituting u = w + uq into the right-hand side of fl2.12p . and by (12. lip , we 
get 



r / \w\^dx= / tr{a{w)Vw)dx - / tr {a d uo) dx + E{i). (2.13) 
Jn Jn Jn 

And the first identity (12.11) follows from (12. 6p . 

Similarly, by substituting uq = u — w into the right-hand side of (12.121) we will 
obtain ([22]) • □ 

3 The testing boundary data 

In this section we construct the boundary data to be input into E for detecting 
the location of dD. For this purpose, we first introduce the CGO solutions with 
complex polynomial phases. 

3.1 CGO solutions with complex polynomial phases 

We are to construct CGO solutions with complex polynomial phases to 

V ■ (To(v) + /cV = (inM^). (3.1) 

Suppose that v G C°^(M^)^ satisfies the above eqaution. By Helmholtz decomposi- 
tion, we can write 

V = + V^ip 

for some smooth scalar functions ip and ip, where V'^ip •= {~d2'4^, diip)^ (and here 
we also regard Vf as a column vector). Then ip and ip satisfy 

V((Ao + 2/io)Ayp + fcV) + V^(;UoAV^ + k^i;) = 0. 
/ .2 \i/2 /l.2\V2 

Let ki = ( Ao+2/Jo ) ~ j ■ -^^^^ above equation it's easy to see 

that conversely for any (p and ip G C°°(M^) satisfying 



A + kjip = 
Aip + k^^/j = 



2„, _n (3-2) 
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V = V<^ + V^V is a solution to (13. ip . Moreover, if ip and ip are CGO solutions to 
( IX^ . then V is a CGO solution to dXT]) . 

It is not difficult to construct CGO solutions to (13.21) by using the Vekua trans- 
form, which transforms a harmonic function to a solution to a Helmholtz equation. 
Precisely, for any real constant w, the Vekua transform T^^ associated with u is 
defined as follows: 

T^(u)(x) = n(x) - ^ ^^(^x)^ { Jo(w|x|v/r^)} rft 

for a function m, where Jq is the zero order Bessel function of the ffist kind. If u is 
a harmonic function, then T^{u) satisfies 

A (TM) + ^\TM) = 0. 

This formula is derived by I. N. Vekua. One can refer to [20] for details and other 
related results. 

In the following we adopt the same idea as in [12] and [13] to construct CGO 
solutions with complex polynomial phases. Given iV G N and /5 G C with = 1, 
let p = Pn,i3 be the function on defined by 

p(x) = /3(xi + zx2)^, (3.3) 

which, by regarding as the complex plane C, is a complex polynomial. Then we 
define 

r = Tm,/, := {r(cos^, sin^) : r > 0, |^ - ^o| < ^} , (3.4) 

the open cone with axis 9 = 9q and open angle vr/A^, where is such that (3 = e"*^^". 
Let r = rjv,/3 := -Re{pAr_^}. Note that in T we have 

r(x) =r^ cosN{e -Go) > 0, 

where x = r(cos 6, sin 6). 

Now for any constant h > 0, eh is harmonic (since it is holomorphic by regarding 
as C), and hence 

= Lff, := Tk^{eh) and ip = ip,, ■= Tk^ie^) 
satisfy (13. 2p . Moreover, and iph are CGO solutions. In fact, we can write 

iPh{^) = e'^(l + Rh,ii^)) and = e'^(l + Rh^^)) (3.5) 

with Rhj (/ = 1, 2) satisfying the following estimates in F: 



\Rk,\ < hf"" 



4r(x 



- 4r X 2r X ' 



(3.6) 



8 



These estimates are established in Lemma 2.1]. In this study we will also need 
estimates of the second derivatives of Rh,i, which are not hard to derive in the same 
manner as the derivation of fl3.6p given in [13] . Actually, by repeatedly applying the 
following well-known recurrence formulae 

where Ji is the Bessel function of the first kind of order 1, and using the basic 
estimates 

\Mt)\<^-, |Jo(t)l<i, vt>o, 

the verification of the following estimates are direct (although somewhat lengthy): 



dxidxj 



< 



2N 



h 

+ 

+ h 



4r^(x) 
kfN{N-l)\x^^ 



4rjv(x) 



+ 



2rjv(x) 



(3.7) 



+ 



kfdij 



in r, for 1 < l,i,j < 2, where 6ij is the Kronecker delta. 

Let diam{Q) denote the diameter of Q. By 03.61) and (13.71) there exists a constant 
Cr = Cr{\o, fiQ, k, N, /3,dia'm{fl)) > such that for any < /i < 1, 1 < /, i, j < 2 
and X G r n 



\Rh.iU)\ < h- 



C 



R 



OEM 



dxj 



< 



< 



TAr(xj 

Cr _ 
^iv(x)' 
1 Cn 



{31 



/irAr(x) 



Now V = v/i := Viph + 'V'^iph is a CGO solution to (13. ip . can be written down 
explicitly as follows: 

Qh,i{^) 

Q/t,2(x) 



P(x) 

e h 



where 



<5h,i(x) 



l ^P(x) ^^ I R 

h dx2 



dRhAM 



x)) + 



dxi 
dRhA^ 



dxi 



(3.9) 
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and 



1 dpjx.) 9i?ft,i(x 
1 9p(x 



+ 



/i dxi 



[1 + RkA^)) + 



dX2 



dxi 



(3.10) 



Thus for < h < 1 and i = 1,2, from (13.81) we have the following estimates for Qh,i 
in r n fi: 



\Qh,i{^)\ < 



N-l 



h 



1 + h 



C 



R 



2Cr 

t{x)J ' r(x) 



+ 



< —tl J tL 

- h r(x)^ 



(3.11) 



and for j = 1,2 



dxi 



< 



N-2 



h 



(iv+ixD^ + iv' 



Cr 



1 + 



r X 



+ 



R 



r X 



2C 



R 



r X 



(3.12) 



where Cr = Cr{Xo, /iq, k, N, (3, diamiVt)) > is a constant. 



3.2 The testing boundary data 

Note that from the discussion above the CGO solutions \h are controllable in F fl fi. 
In the following we go on to follow the idea in [T^] and [13] to modify into a 
family of functions localized in F. 
For t > 0, let 

:= {x G F : r(x) = i}, (3.13) 

the level curve of r in F at j. In fact, any level curve of r = r^r^^ has branches, 
and the cone F = F^r,/? just contains one branch with the two edges of F being the 
asymptotes of that branch. Also note that when t is larger, the curve it is closer to 
the origin. (We refer to Figure 2 in [13] for an illustration.) Then for d > let 

Td = TN,M-= U (3.14) 

0<t<d 

Note that for di > d2 > we have F^j C F^^. 

In the following we fix an e > and a compact interval J C (0, oo). Let {4>d}deJ 
be a family of smooth cut-off functions such that 
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(i) < 0d(x) < 1, 

(ii) 0d(x) = 1 (resp. 0) for x G Td+e (resp. x G \ Td+2e), and 

(iii) for some > 0, we have \d"(j)d{y^)\ < for each multiindex a with \a.\ < 2, 
for each x G and for each d E J. 

The existence of such family {(pd} is obvious and we omit a precise construction. 
Now let 

Pd,Hi^) := 0d(x)e-7^v^ G ^"^(M^)^. (3.15) 

It is the traces of these pd,h on dfl that will be the testing data to be input into E. In 
fact, we will see that the behavior of E{pd,h\dn) as /i — )■ O"*" tells whether Td intersects 
D or not. Now note that although pd,h is controllable from the discussion above, it 
is no longer a solution to fl3.ip . However, to get information from E{pd^h\dn) we will 
need estimates related to the solution of (13. ip with boundary condition pd,h\dn- But 
indeed for small h controllability of pd,h gives controllability of the true solution of 
(13.11) with the same boundary condition. We explain this precisely in the following. 
Let uo d,h satisfy 



V ■ cro(uo,d,h) + k'^Uo^d,h = inQ 
uo,d,h = Pd,h\dn on dQ. 



(3.16) 



And let 

Wft = Pd,h - uo,d,ft, (3.17) 

then Wh satisfies 



V ■ ao(w;,) + k'^Wh = V ■ (To(pd,ft) + k'^pd,h in fl 
Wh = on dQ. 



(3.18) 



Let 

gh = ^ ■ (roiPd,h) + k'^Pd,h: (3.19) 
then we have the following lemma. 

Lemma 3.1. There exists positive constants Ci and C (depending on Q, Xq, fiQ, k) 
such that for < h < 1 and d E J 

C 111 

In particular, there exists < Hq < 1 such that for < h < ho and d E J , there is a 
positive constant C = C'{Q, Ao,/io, k) such that 
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Proof. Note that in this paper for any two vectors a and b, we define a ® b to be 
the matrix whose ij-th entry is Qibj. 

That ||w/j||//i(n)2 < Ci||gfe||i2(Q)2 for some Ci is classicaL So we need only to 
estimate ||g/i||L2(s^)2. 

Since Pd,h(x) = 0rf(x)e~Mv,„ we have 



gh = e-^|AoV(V ■ {(j)dVh)) + /ioV ■ (V(0dV;,) + (V(0rfvO)^) + k^(l)d^h 
= e-^|Ao [V(V0rf ■ V;,) + V0d(V ■ v^)] 
+ /ioV ■ [vh V(pd + V0d (g) Vh] 

+ /^o(Vv;, + (Vv;,)^)V0rf 

Because V ■ o"o(v/i) + k^^h = and Vcpd = outside Td+2e \ ^d+e, we have 

WShh^n)^ < C'se"^||vft||j^i((r^^2^\rrf+,)nf7)2 (3.20) 
for some positive constant Cg = Cg{\Q, hq,C^). 

By dsm), for X G (r,+2. \ r,+,) n fi, 



|v,(x)|=e"^^|Q,,i(x)|2 + |g,,2(x)|2 

for some positive constant C'^ = C^(Ao, /io, k, diarriQ). Hence we have the following 
estimate: 



h 

Similarly by (13. lip and fl3.12p we have 



2t(x) 

e 

(rd+2,\rd+,)nc 



V-'(rd+2Ard+a:)nn 
Since 

/ e'-^dx <\(Td+2e\Td+s)nn\eT^^, 

we have 

1 C 1 1 1 

||g/.||L2(f7)2 < C'9e'^l|v/,||//i((rd+2Ard+.)nc)2 < ;^e"^^^~^\ (3.21) 

where C depends only on Aq, Ho, k and Q. 

□ 
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4 The main theorem for the reconstruction of un- 
known inclusions 

We now come to considering our inverse problem of reconstructing D. For the main 
theorem we make the following three assumptions (in addition to those already made 
in the introduction) throughout this section. 

1. We assume V ■ cxo and V ■ a satisfy the strong convexity condition (but not 
only the strong elliptic condition (11. 3p ): 

Ao + yUo > 0, /io > 0; 
A + /i > 0, yu > on fi. 

Thus, in particular, Lemma 12.21 applies. 

2. (Ad + /iD)/iD > on D. 

3. For any y G dD, there exists a ball -Br(y) such that one of the following jump 
conditions holds: 

(i) ^d(x) > r, Az5(x) + /iD(x) > 0, Vx G 5,(y) n 

(ii) ^Z3(x) < -r, Az5(x) + < 0, Vx G 5,(y) n ^ ' ^ 

Now assume the origin is outside As in section 3, in the following we fix 
an G N, a /3 G C with = 1, an e > 0, and a compact interval J C (0, oo). And 
recall the definition of p, F, F^, and pd,/, in ([33D, dM]), (13131) . (l3AiD . and (l37[5|l 
respectively. Also recall that we use r to denote Re{p). Let 



sup r(x), if DnF ^ 
xeDnr 

, if D n F = 



Note that -D fl F 7^ if and only if > 0, and in this situation ii/s^ is a curve just 
touching dD, i.e. iys, nD = iys, DdD^^. 

For notational simplicity let id,h '■= Pd,h\dn- Recall that UQ^d,h satisfies 

V • (To(uo,d,ft) + k^uo^d,h = inn 
uoA,h = U,h on dQ. 

Similarly let be the solution when the inclusion D exists: 

V ■ (T{ud,h) + k'^Ud,h = 0, inQ 
^d,h = U,h on dn. 



general, for a = (ai, 02)"^ a point outside fj, we should use p = f3{{xi — ai) + i{x2 — 02))^, 
and similar modifications of F, Td, etc., and there is a similar result as Theorem 14.11 However as 
we can always set the coordinates so that ^ f2 in practice, such consideration is not needed. 
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Now let Wrfjj = Urf /j — UQ^d,h- We have the following two inequalities from Lemma 



E{id,h) < / (Ad + /iD)|V ■ UQ^d,h\'^dx 
Jd 



2 / /iD 

'd 



(4.2) 



e(uo,d,h) - 2(^ ■ ^OAh)h 
'd A + /i 



+ 2 



D 



e(uo,(i,h) - 2(^ ■ Uo,(i,/i)/2 



(4.3) 



They are the key to the following main theorem of this paper. 

Theorem 4.1. For d E J and h > small enough, the following conclusions hold: 

(A) IfDnTd = (lS, then 

for some C > independent of h, where Sd = max(^^, s=k) < ^. 

(B) If D r\T d ^ ^ and D has continuous boundary, then there exists a constant 5, 
< 5 < — such that 

md,h)\>Ch-^ei^'--^-'^ 
for some C > independent of h. 
(B') If DnVd^% andD has C"-" boundary for\<a<l, then 

|i?(W)| >C/i-3+ie^(^*-H) 
for some C > independent of h. 



Before going into the proof of Theorem 14. H we give two remarks. 
Remark 4.1. 

1. From Theorem 14.11 we have the following conclusions. In (A), since Sd < ^, 
\E{id,h)\ tends to zero as h tends to zero. On the other hand, in (B) and in 
(B'), since s* > ^, |_E'(fd,/i)| tends to infinity as h tends to zero. In particular, 
from (A) and (B), we have 



(4.4) 



= inf <^ - : lim \E{fd,h) 
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Hence, although we don't know the hmiting behavior of E{iii h) when just 
touches dD, we can reconstruct dD in principle. (Of course, due to the ge- 
ometric nature of F^, in fact only "detectable" points can be reconstructed. 
An explanation of this point can be found in [18, Corollary 5.4]. Also see [19] 
or [13] for a reconstruction algorithm, which is easily modified to be suited 
for our case. We omit such discussions in this paper.) From this point of 
view, almost no regularity assumption on dD is essential in the reconstruc- 
tion. Nevertheless, for a complete characterization of the limiting behavior of 
E{fd,h), we include (B') in our theorem, while for this purpose more regularity 
assumption has to be made. 

2. We will use f l4.2p and fl4.3p to prove Theorem 14.11 Roughly speaking we have 
better knowledge of uo^d,h than Wd,h, and the crucial step is to give an ap- 
propriate control of ||wd^/i||x,2(Q) in terms of UQ^d,h- For this purpose, in the 
corresponding parts of early researches, e.g. [6l HH |T3], some technical as- 
sumptions (precisely, positivity of the relative curvature and finiteness of the 
number of touching points of ii/s^ (or say Fi/^J and dD) have to be made. In 
particular dD is usually assumed to be C^. (In order to apply CGO solutions 
with complex polynomial phases, even more technicalities are involved. For 
example, in [131 Lemma 3.7], the authors proposed an estimate which is based 
on a rather technical result in [ID]-) In PLSJ, Sini and Yoshida came up with a 
totally different method to control ||wrf /j 11^,2 (j^) (while they did not adopt CGO 
solutions with complex polynomial phases). Precisely, they proposed (in our 
terminology) 

||Wd,/i||L2(Q) < C||uo,d,/i||iyi.P(D) (4.5) 

for some p < 2, which was proved by using an regularity estimate of 
Meyers and the Friedrichs' inequality. In this way the technical assumptions 
on the touching point are no more needed and dD can be assumed to be 
only Lipschitz. Inspired by this result, we tried to adopt their idea in our 
situation. We find it's interesting that, with more careful treatment, we find 
the boundaries of the inclusions can in fact be assumed to be only continuous. 
Moreover, we find in the case of F^ just touching dD, the regularity assumption 
on dD can be reduced to be C°'" for any a G (|, 1]- 

To save notation, in the remaining of this paper we will freely use C to denote 
a constant, which may represent different values at different places. 

The following lemma is just (14. 5p . we give the proof here for the sake of com- 
pleteness. 

Lemma 4.2. There exist constants C > and 1 < < 2 such that for qq < q < 2, 

||w||L2(f^) < C|| Vuo||l'?(D), 

whenever u and Uq G H^(Qy satisfy (11.11) and (II. 2p respectively, u and Uq have the 
same traces on dfl, and w = u — uq. 
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Proof. Let q be the element in HqIQ)"^ satisfying 

V ■ (o-(q)) + fc^q = w mQ. 
Then, by taking inner product with w and integration by parts, we have 

/ |w|^(ix = — / tr(cr(q) Vw)(ix + A;^ / q^-wdx 
Jq Jq Jn 

= — tr{a{w)Vq)dx + k"^ / q-wdx. (4.6) 
Jn Jn 

On the other hand, note that 

V ■ cr(w) + k'^w = —V ■ cr /)(uo), 
which, by taking inner product with q and integration by parts, gives 

— / tr{a{w)\/q)dx + A;^ / w ■ qdx = / tr(a£)(uo)Vq)(ia;. (4.7) 



From 04.61] and fl4.7p we get 

/ |wp(ix = / tr{a£){uQ)\/q)dx. 
Jn Jn 

Then by Holder's inequality we have for any 1 < p < oo 

/ |wprfx < ||(T£,(uo)||L9(D)||Vq||LP(n), (4.8) 
Jn 

where q is the conjugate exponent of p. 

Now let Q = q. By definition of q we have 

V ■ (a(Q)) = w - A;2q in 
Q = onQ. 

Then by [121 Theorem 1], there exist Pq> 2 such that for each 2 < p < po, 

||Vq||LP(n) = ||VQ||LP(n) < C{||q||L2(j^) + ||w||L2(f^)} (4.9) 

for some C = C{k, A,/i) > 0. Note that also by definition of q, we have ||q||i2(f^) < 
C||w||l2(q) for some C = C{k,\,fi) > (see e.g. [21 Section 6.2, Theorem 6]). So 
from (14. 9 P we have 

l|Vq|Up(n) < C||w|U2(f,) (4.10) 
for some C = C{k, A,/i) > 0. Combining (14.81) and (14.101) . we have 

||w||^2(Q) < C||Vuo||L9(D)2||w||L2(f^) 

for some C = C{k, X, fi) > and 2 < p < pq, and therefore 

||w||i2(Q) < C|| Vuo||l<7(d) 

for some C = C{k, X, fi) > and qo < Q ^ 2, where 1 < go < 2 is the conjugate 
exponent of po- D 
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Remark 4.2. Remember that we assume has a smooth boundary for simphcity. In 
fact, it is so assumed only to allow direct application of the estimates in [12]. In 
other words, the regularity condition on dfl really required is just that guarantees 
the validity of f lOj) . 

To make the proof of Theorem 14. II more concise, some computational results are 
collected in the following lemma. 

Lemma 4.3. For d & J, we have the following conclusions. 

(i) There exsists a constant C > such that for g > and Q < h < 1, we have 

\\yVdA\L^{D)<Ch-^\ I e^(^W-^)da;) . (4.11) 

In particular, since s* > r(x) for x G -D fl Td+2e, we have 

||Vp,,,||^,(^) < Ch-'e-h^'*--^^ 
for some C > independent of h. 

(ii) There exist positive constants c and C such that, for < /i < 1 and for any 
open set U with D fl Td+e H U (/), we have 

1 ' 
e{pd,h) - -(V ■ Pd,h)l2 

> {ch-^ - Ch-^) [ e^("(")-3)rfx. 

JDnTa+^nU 

(iii) There exist constants C > and Qq < 2 such that for each qo < q < 2 and 
< h 1, we have 

\\^d,h\\Un) < Ce^(^-^) + Ch-^ ( [ e^("W-^)cia;) . (4.13) 

Proof. (i) Remember that 

P<i,h = {pIh.pIhT = 0.ei(^W-^)(Q,,i,Q,,2)^, 

where Qh,i and Qh,2 are defined in (13. 9p and (13.101) respectively. Then by 
definition of (pd (in page[TT]), for x G -D \ Td+2£, we have prf^/i(x) = 0. On the 
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other hand, by (I3.1ip and f l3.12p . we have for x G fl Tii+2e and < h < 1, 



dxi 



dxi 



(4.14) 



+0d X 



h dxi 



dxi 



for some positive constant C independent of h. Since s* = sup r(x) and 

xGDnr 

iVpd,/,!" = dVp^.h!')"/', we have for < /i < 1 



llVpd, 



h\\Li{D) 



'Dnrd+2, 

for some positive constant C independent of h 



(ii) We can compute ^ directly by and fIXTU]) for x G D n r^+ae: 



dxi 



and 



dxi 



_J_ p(x) 

e 



e 



_ j_ 

e hdQ h 



dx 



h dxi 
dp{x.) 9p(x)\ dp{x.) 



\ dxi 



dxi 



dxi 



dxi 



Qh,2 + TVd^ Qh,2 + (pd- 



dxi ^"'^ ' h^"" dxi 
j_ pW / 1 / f^p(x) (9p(x)\ (9p(x) 



\ dxi 



dxo 



dxi 



dxi 



where 



1 dp f dp 



dxi \dxi 



-RhA — 



dp 
dxo 



+ 1^ 
h dxi 



(i + fi 



dp 
dxi 

dp fdRh,i dR, 



h,2 



dp 

dxo 



{l + R 



h,2) 



h,2 



+ 0d 



dxi \ dxi 
d^p 



dxo 



+ 



dp dRh^i dp dR, 



h,2 



dxi dxi dx2 dxi 



(1 + Rh,i) - 



dxidxi 
d'^Rh,! _ d'^Rh,2 
dxidxi dxidx2 
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d^p 
dxidx2 



{l + R 



h,2 



+ 



dRhi dRh2 



dxi \ dxi 



dxo 



1 dp f dp 



dxi \dx2 

d 



RhA + 



dp 
dxi 



R. 



h,2 



1 

+ - 



h dxi 



8X2 

dp f dRh,i 
dxo 



+ 



dxi 
d'^p 



+ 



dxi 

dRh,2 
dxi 



^ dp dRh,i _^ dp dRh,2 
dx2 dxi dxi dxi 



(1 + Rh^) + 



dxidx2 

d'^Rh,! _^ d'^Rh,2 



+ 



d^p 
dxidxi 
d(j)d 



{l + R^ 



h,2) 



dxidx2 dxidxi / dxi \ dx 
By dSSD, for any X G D n Vd+2e and < /i < 1, 

l4-i(x)|,|j;-.(x)|<c/.-i 

for some positive constant C independent of h. 
Then we have for x G -D fl Td+2e and < /i < 1 



e iT>d,h) - ■ Vd,h)l2 



dRh,i _^ dRh,2 



dxi 



> 



> 



d,h 



dp f dp dp 



dxi dx2 



> ei^^--^\c<Plh-^ - 2Ch 



dxi \dxi dx2 

2 

h 

~2 



dp 

dx2 \dx2 ' dxi 



dp ^ dp 



for some positive constants c, C independent of h. Then (ii) of this lemma is 
vahd. 



By Lemma [4. 2[ there exist constants C > and 1 < go < 2 such that 

||Wd,/i||L2(!^) < C|| Vuo,d,/i||L9(D) 

for each qo < q < 2. Therefore replacing UQ ^ h by Pd,h — and applying 
Holder's inequality, we have 

||Wd,/,||L2(f^) < C {\\VWh\\Li(D) + ||Vpd,ft|U<7(Z))} 
<C{\\VWh\\LHD) + ||Vpd,J|M(^)} 
< C {\\VWh\\m{D) + ||Vpd,/,||L<,(D)} • 

Then by Lemma lO and fICTD . KT^ follows. 

□ 



19 



Now we give the proof of the main theorem. 
Proof of Theorem \4.1\ (A) By (14.31) . we have 

Jd a + zh 

2 

+ 2 / e(uo,,,„)-^(V-uo,,,,)/2 
Jd ^ 



dx 



Therefore, together with (14.21) . we have 



for some positive constant C independent of h. Therefore from Lemma 13.11 
and Lemma [4.31 we have, by choosing q = 2, the following estimate: 

|^(fd,/^)| <C|||VW,|| 

L2(n) + llVpd ,/i|Il2(£)) + ||Wa!,/i|||2(Q)| 

<c|et(^-^) + /i-4ef(--^)|. 
Therefore for < h < 1, 

where Sd = max(^, s*). Moreover we notice that Z) fl = implies < |, 
and the conslusion (A) follows. 

(B) We first consider case (i) of f l4.1l) and prove the conclusion (B) from fl4.3p . 

Suppose D n Frf 7^ 0, then s^, > ^ > ^ since D is open. Therefore for any 
y G dD n Td+e, each neighborhood Uy of y satisfies D fl Td+e H f/y 7^ 0. By 
the assumption (i) of (14. ip . for each y G there exists Vy such that 

/iD(x) > Ty, Xd + I^d> 0, Vx G fi,^ n D. (4.15) 

Set K := dD n {t = s^} = dD n £1/,.. It's easy to see that K ^ (/). Since 
K is compact and is contained in U^^gx-Bry (y), there exsits G N such that 
K C Uf=i5,^(yj), where Ty^ is abbreviated to rj. Let Dr — D \ U^LiBr^iy 
then it is easy to see that there exists 6' > such that 



r(x) <s^-6' in 



R- 
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Therefore for qq < q < 2 we have 

N 

e 



JDnr^.o. Jdr ,_i Jb, 



'Br, (y,)nDnrd+2^ 



for some y* G {yj}^i and r* G {rj}^^ such that 



e 



d)(ia; = max ( / e^^'^ d^dx 



Moreover, we can compute more finely that 



e^^^ d^dx 

e^^'^~^^dx+ / e^^^~^^rfx 
Br, (y* )nDnrd+, J Br, {y, )nDn{ra+2eVd+e ) 



< / e'^*^'^ d^dx + Ce'^'''d+e d\ 

J Br, (y.)nDnrd+e 
Therefore by combining the above inequahties, we have 

ih^'^^i^dx 



[ e^^^'-d^dx <C [ 



'Dpra+2e J Br,{y,)nDnra+e 

+ Ce^.(iih--d^ + CeTA^'--d-''\ 

Set 

9;* 



A.,, := I e^^^-^^d. 



X. 

Br,{y,)nDnra+e 



Now we come back to (14.31) , from Lemma 13.11 we have for < /i ^ 1 

flQfiD 



D 
2| 



1 ' 

4Pd,h) - -(V ■ Pd,h)l2 



dx-\\wh\\Hi(n) 



> C 



Id 



1 ' 

g h d+e d I 



dx 



(4.16) 
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In the following we estimate each term separately. 
First, by Lemma [4.31 we can compute 



12 



> c- 



A2,.,hich-'' - Ch^) 



(4.17) 



for qo < q < 2. 

Now we need to compute Ag^^: h carefully. For y^, G ii/s, H dD, we consider 
the following change of coordinates as in [13j. First, let T be the composition 
of the following two rigid motions: i) translate y,,, to the origin, and ii) rotate 
so that the unit inward normal of T{Ti/s^) at the origin is the vector (0, 1)"^. 
Then set z = (2;i(x), 2;2(x))-^ = T(x) and ^ = (^i(z), ^2(2))"^ = 2(z), where 



^ z 



Zl 

r(T"^z) - 



Then S o 7" gives a diffeomorphism in a neighborhood Uy^ of y*. Geomet- 
rically, under the transformation H o 7" the point y* becomes the origin of the 
new frame (^1,^2)^5 (^i-axis coincides with the curve ii/s^, and the positive 
direction of ^2-axis coincides with the unit inward normal of T(ri/sJ at y^,. 

We do the above change of coordinates, then we have S o T(y*) = and 



\J~oTiBr, (y.)n-Dnrd+,) 

SoT{Br,{y*)nDnTa+,) 



(4.18) 



Since dD is continuous, S o T{dD) is also continuous and is able to be 
parametrized by a continuous function near ^ = under a suitable rotation. 
So, we consider a rotation T with T(^) = ^ = {^1,^2)^ such that T(HoT(9D)) 
can be parametrized by /*(^i) near ^ = with /*(0) = 0. 

Actually, we can choose T such that 

^2 = (sin^)^! + (cos^)f2 with |^| < ^, 

because S o T{D) C {^2 < 0} and D is open. Let a = sin 6* and h = cos 9, 
then b > 0. Without loss of generality, we assume T(S o T{dD)) can be 
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parametrized by in < diam{T{^ o T{BrXY*) H D fl Vj_^^))). Set 

t/ = f(Sor(5,,(y,)nDnr,+,)). 

Here we note that /* is continuous in U since we assume dD has continuous 
boundary, and U C {a(i + < 0} since S o T(-D) C {^2 < 0}. 

Now it's easy to see that there exist positive constants 6i,62,S[, 62 independent 
of h with 62 < S2 such that 

J-5[ J -52 

JsoT{BrAy*)nDnrd+e) Ju 





'-5[ J -52 

Since t/ C {a^i + b^2 < 0}, 

S2<mi)<-^ii in ^. 
Therefore, we can compute directly and obtain that 

< / e^^M^ (4.19) 

= / / e^(<^+'«^)rf|+(5; + 5i)(5^-52). 
In order to make the computation clear, we set 

Bg,,,H= / / e*(<^+^«^)d|. (4.20) 

J-S[ J -52 

By combining fl4.17p . f l4.18p and f l4.19p . we have 
I|wd,^lli2(n) 

. ^ e^^"-^^i?2,.,/.(c-C/^^) .4 

— ^ 2/ In, ^xn; 2/1 In o , 1 V • / 



c 



^2,.,fe(c-C/i^) 

(5g,.,^)2/5 + C + e^(^-^-) + et(-^') ■ 
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Now we compute more carefully. We note that since is continuous 

near ^ = 0, for all < 5 < min((52, s*— 6'), there exists < 6'1 < min(5i, 6[) 
such that 

= \Mii) - /*(0)l < S, V6 e i-5'l,5'l). 

Therefore, 



Bn,*,h > 



> 



Mil) 



-52 



-5" Qb 



> — e h 
qb 



Then for < /i <^ 1, we obtain the following estimate 

h _q - r 

Bq,*,h > C—e h 



qb 



(4.22) 



for all < 5 < min(52, s^—-^, 6') and for some C independent of h. Moreover, 
we observe that for < /i <^ 1 



B. 



and 



2,*,h 



e h° 



2,*,h 



Then (14.21 1) becomes the following estimate 



> c 



B' 



2,*,h 



for go < g < 2 and < /i < 1. 

Actually, we can directly compute and use the Holder inequality to obtain that 
for qo < q < 2, 



-5' Qb 



2/q 



qb) 



2/q 



Si 



2/q 



2/q pSi 



2/g 
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Since 82 + /*(^i) > for ^1 G [—5^, ^i], we have 

< e-^(^^+^*(«^» < 1 

and therefore for g < 2 



1 _ g-f (<52+/*(Sl))y^'' < 1 _ g-f (^2+/*(Sl))_ 



Hence we obtain for qo < q <2 

< C (J^y f' e*°«'ei'-'-'«'' (1 - (if 1 



g6/ \h 




h 



2/q 



2b 



and then the most difficult part of the proof of Theorem 14.11 can be concluded 
that for < /i ^ 1 and qo < q < 2 



H 112 >Ch 1, 

for some constant C independent of h. 

Back to dnS]), by we have for < < 1 and qo < q < 2 



(4.23) 



E{fd,h) > C 



D 



^{Vd,h) - ■ Vd,h)h 



dx 



X 1 



2(^^ 1) 



(4.24) 



By direct computation we have 



therefore 



/D^h(P'^,^)-|(V-P.,.K2rcia; 



2 1 „ ^ 



0(1). 



(4.25) 
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Hence by using Lemma 14.31 and by computing directly from fl4.24p and fl4.25p , 
we have for < /i <^ 1 

|^(fd,h) I > Ch-'A2,.,H > Ch-'e-^^''--^^B2,.,H. (4.26) 

Therefore by fl4.22p . for all < 5 < min(52, s*-^, S') and for < /i < 1 we 
have 

\E{Uh)\>Ch"'e-^^''--^-'\ (4.27) 

for some constant C independent of h. Choose 6 such that 6 < s* — ^, then 
the proof of (B) is complete. 

For case (ii) of (14.11) . instead of using (14. 3p . we shall consider the negative of 
(1121): 



-E{id,h)> / -(Ad + /iz))|V ■ uo.d,/,prfx 
Jd 



2 / fiD 



D 



1 ' 



dx — k'^\W(i hl'^dx. 



n 



And a similar argument will also give (I4.27p . 

As in (B), we will only prove case (i) of (14. ip by (14. Sp . and case (ii) of (14.11) can 
be treated similarly by using the negative of (14. 2p . Suppose that DnTa ^ and 
D has boundary. Since DnTa ^ 0, > ^ > and K = dDnii/s, ^ 0. 

In fact, we have proved in (B) that s^: — -^ > and continuity of dD ensure 
(I4.23P holds. So (I4.23p also holds under this assumption of (B') and therefore 
flCTj) also holds. 

However, since D has C°'" boundary, we have the better estimate than (14.221) . 
Without loss of generality, we assume /*(^i) is C°'" for G [— Then 
there exists a positive constant L such that for G [—6 — 

\Mii)\ = \M^i)-fM\<L\^ir- 
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Therefore we can compute directly as follows: 

J-S[ J -52 
J-5[ \J-S2 







2a ; 



2b 

Without loss of generality, we assume that < Si < 1. Since < a < 1, we 
have 

Jo Jo 

Jo 

Then by computing directly, we have 



Since 

4 



and 



Jo Jo 



e-^'^^'^fi ^ / e-^'^^'^li < OO as /i ^ 0+, 



we have 

52,*,h > C/i^+i for < /i< 1. 

Then by fOB]) 

for each < /i <^ 1. If a > |, then even when s,, = ^, \E{fj_^h) \ tends to infinity 
as h tends to zero. 

□ 



27 



Acknowledgements 



I would like to deeply thank Professor Jenn-Nan Wang for the suggestions of this 
problem and for many helpful discussions. I would also like to thank Professor Gen 
Nakamura for informing me about the work of Sini and Yoshida. Thank Lircn Lin 
for many discussions on my English writing. This work is partly supported by the 
National Science Council of the Repubhc of China. 



References 

[1] A. p. Calderon. On an inverse boundary value problem. In Seminar on Nu- 
merical Analysis and its Applications to Continuum Physics (Rio de Janeiro, 
1980), pages 65-73. Soc. Brasil. Mat., Rio de Janeiro, 1980. 

[2] Lawrence C. Evans. Partial differential equations, volume 19 of Graduate Stud- 
ies in Mathematics. American Mathematical Society, Providence, RI, 1998. 

[3] T. Ide, H. Isozaki, S. Nakata, S. Siltanen, and G. Uhlmann. Probing for 
electrical inclusions with complex spherical waves. Comm. Pure Appl. Math., 
60(10):1415-1442, 2007. 

[4] M. Ikehata. Reconstruction of the shape of the inclusion by boundary measure- 
ments. Comm. Partial Differential Equations, 23(7-8) :1459-1474, 1998. 

[5] M. Ikehata. Enclosing a polygonal cavity in a two-dimensional bounded domain 
from Cauchy data. Inverse Problems, 15(5):1231-1241, 1999. 

[6] M. Ikehata. How to draw a picture of an unknown inclusion from boundary 
measurements. Two mathematical inversion algorithms. J. Inverse Ill-Posed 
Probl, 7(3):255-271, 1999. 

[7] M. Ikehata. The enclosure method and its applications. In Analytic extension 
formulas and their applications (Fukuoka, 1999/Kyoto, 2000), volume 9 of Int. 
Soc. Anal. Appl. Comput., pages 87-103. Kluwer Acad. Publ., Dordrecht, 2001. 

[8] M. Ikehata. An inverse transmission scattering problem and the enclosure 
method. Computing, 75(2-3):133-156, 2005. 

[9] M. Ikehata and H. Itou. An inverse problem for a linear crack in an anisotropic 
elastic body and the enclosure method. Inverse Problems, 24(2):025005, 21, 
2008. 

[10] Y. Li and L. Nircnberg. Estimates for elliptic systems from composite material. 
Comm. Pure Appl. Math., 56(7):892-925, 2003. Dedicated to the memory of 
Jiirgen K. Moser. 



28 



[11] W. McLean. Strongly elliptic systems and boundary integral equations. Cam- 
bridge University Press, Cambridge, 2000. 

[12] N. G. Meyers. An L^-estimate for the gradient of solutions of second order 
elliptic divergence equations. Ann. Scuola Norm. Sup. Pisa (3), 17:189-206, 
1963. 

[13] Sei Nagayasu, Gunther Uhlmann, and Jenn-Nan Wang. Reconstruction of pen- 
etrable obstacles in acoustic scattering. SIAM Journal on Mathematical Anal- 
ysis, 43(1):189-211, 2011. 

[14] G. Nakamura and K. Yoshida. Identification of a non-convex obstacle for acous- 
tical scattering. J. Inverse Ill-Posed Probi, 15(6):611-624, 2007. 

[15] M. Salo and J.N. Wang. Complex spherical waves and inverse problems in 
unbounded domains. Inverse Problems, 22(6):2299-2309, 2006. 

[16] M. Sini and K. Yoshida. On the reconstruction of interfaces using ego solutions 
for the acoustic case. Preprint. 

[17] G. Uhlmann and J.N. Wang. Complex spherical waves for the elasticity system 
and probing of inclusions. SIAM J. Math. Anal, 38(6): 1967-1980 (electronic), 
2007. 

[18] G. Uhlmann and J.N. Wang. Reconstructing discontinuities using complex 
geometrical optics solutions. SIAM J. Appl. Math., 68(4):1026-1044, 2008. 

[19] G. Uhlmann, J.N. Wang, and C.T. Wu. Reconstruction of inclusions in an 
elastic body. J. Math. Pures Appl. (9), 91(6):569-582, 2009. 

[20] I. N. Vekua. New methods for solving elliptic equations. Translated from the 
Russian by D. E. Brown. Translation edited by A. B. Tayler. North-Holland 
Publishing Co., Amsterdam, 1967. North-Holland Series in Apphed Mathe- 
matics and Mechanics, Vol. 1. 



29 



